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Filling-control metal-insulator transition on the two-dimensional Hubbard model is in- 
vestigated by using the correlator projection method, which takes into account momentum 
dependence of the free energy beyond the dynamical mean-field theory. The phase dia- 



H 

gram of metals and Mott insulators is analyzed. Lifshitz transitions occur simultaneously 
with metal-insulator transitions at large Coulomb repulsion. On the other hand, they are 
separated each other for lower Coulomb repulsion, where the phase sandwiched by the Lif- 
shitz and metal-insulator transitions appears to show violation of the Luttinger sum rule. 
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Through the metal-insulator transition, quasiparticles retain nonzero renormalization factor 

(\ ' and finite quasi-particle weight in the both sides of the transition. This supports that the 

metal-insulator transition is caused not by the vanishing renormalization factor but by the 

relative shift of the Fermi level into the Mott gap away from the quasiparticle band, in sharp 

l/^j ■ contrast with the original dynamical mean-field theory. Charge compressibility diverges at 

zi. ■ the critical end point of the first-order Lifshitz transition at finite temperatures. The origin of 

the divergence is ascribed to singular momentum dependence of the quasiparticle dispersion. 
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*T3 ' !• Introduction 

a: 

Transitions between metals and insulators caused by electron correlation effects have been 
O ■ 

studied extensively. 1 A typical example of the correlation-driven metal-insulator transitions is 



the Mott transition, near which various competing orders including unconventional supercon- 
ductivity in the copper oxides are found. In spite of long history, the nature of Mott transitions 
still remains a challenge of recent active research and a controversy remains. A key issue of the 
Mott transition 1 has already been posed by early studies by Hubbard 2 and Brinkman-Rice, 3 
which have proposed somewhat contradicting scenarios each other. The former describes the 
formation of energy gap and the insulator appears from the shift of the Fermi level into the 
gap, while the latter describes the transition to the insulator by the mass divergence of the 
quasiparticle. 

Recently, the two views are unified in the limit of infinite dimension by the dynamical 
mean- field theory (DMFT). 4,5 DMFT succeeded in describing the two aspects of the Mott 
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transition; formation of an energy gap and vanishing of the Kondo-like resonance at the Fermi 
level. The latter is described as vanishing quasiparticle weight at the Fermi level Z — > 0. 
Although the Mott gap is formed already in the metallic phase before the transition, the 
criticality of the metal- insulator transition is governed by the vanishing renormalization factor 
of the coherent band at the Fermi level. 

Unfortunately, however, DMFT is a theory which is exact only in the infinite-dimension 
limit and its validity in finite dimensions is not clear. In finite dimensions, the self-energy be- 
comes momentum dependent in contrast to DMFT. Results of recent angle-resolved photoe- 
mission spectroscopy (ARPES) measurements 6 ' 7 have revealed momentum dependent aspects 
of metal-insulator transition (MIT) in low-dimensional systems like cupurate superconductors. 
In these systems, momentum space anisotropy of quasiparticles is especially enhanced in low- 
doping samples, which implies that the momentum dependence of the self-energy may play a 
crucial role in MIT and physics around it. Therefore, it is desired to extend DMFT to allow 
the momentum dependence and clarify how it modifies the nature of the transition. 

The central issue is whether the Mott transition in finite dimensions is caused by the 
vanishing of the quasiparticle at the Fermi level or caused by the move of the quasiparticle 
dispersion out of the Fermi level. The answer to this fundamental question may depend on 
the way of treating the momentum dependence of quasiparticles and therefore one has to be 
careful in obtaining reliable answer in finite dimensions. In this paper, we study this issue by 
considering the momentum dependence of the self-energy beyond the dynamical mean-field 
theory and clarify the fate of quasiparticle in realistic finite dimensional systems. 

We reiterate the two possible scenarios of the Mott transition in more detail in the fol- 
lowing: 

(1) Quasiparticle rigid-band picture 

It is based on a certain 'band structure' of quasiparticles. Here, the 'band structure' does 
not necessarily mean the non-interacting one. It is the energy dispersion of quasiparticles, 
which may have gap-like structure formed by electron-electron correlations. MIT is described 
as Fermi level shift out of the top/bottom of the quasiparticle band. Namely, Fermi level 
excitations disappear because the Fermi level moves to the energy range of the gap where there 
are no quasiparticles. Its prototypical example is the Hubbard picture of MIT, 2 although the 
Hubbard approximation itself clearly oversimplifies the real transition. If this route is realized, 
at least the quasiparticle weight closest to the Fermi level remains nonzero all through the 
transition. 

(2) Quasiparticle-weight vanishing picture 

The other type of description is similar to the Brinkman-Rice scenario 3 stating that the weight 
of the quasiparticle disappears at the Fermi level. Quasiparticle poles stay at the Fermi level 
by keeping the Fermi surface until the transition but their weights disappear at the transition. 
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The MIT is purely driven by the weight disappearance. We categorize the DMFT scenario of 
MIT here since it describes the MIT as a formation and a destruction of a resonant excitation 
at the Fermi level. 8 

We note that the real transition may occur in a compromized manner and may have 
the both characters. In other words, in reallity, both of the gap formation and the weight 
renormalization may occur simultaneously. However, there still remains a crucial issue whether 
the quasiparticle weight remains or not at MIT. In its criticality, the two routes are not 
compatible at least when the transition is continuous without the phase separation. 

Another related issue in two-dimensional systems is the relationship of the Mott transition 
to the Lifshitz transition. If the Mott gap is preformed in the metallic phase, the Lifshitz 
transition may trigger a transition from large Fermi surface to small pockets through the 
arc-like structure. However, the first-order Lifshitz transition may preempt this phase with 
the small pockets and cause a simultaneous Lifshitz and metal-insulator transition. The fate 
of the quasiparticle might depend on these different choices of the routes. Such an issue is 
relevant only in finite dimensional systems and has not been explored in detail. 

Recently there have been several attempts of momentum-sensitive analysis beyond DMFT. 
The cellular-DMFT (cDMFT), 9 the dynamical cluster approximation (DC A) 10 and the cluster 
pertubation theory (CPT) 11 have been applied with certain successes. However, the resolution 
of analyzing spatial fluctuations is limited because it is difficult to increase the number of 
clusters in space. This makes it difficult to discuss growth of short-ranged correlations. Here 
instead, we use the correlator projection method (CPM), 12 where spatial fluctuations can be 
considered with a large number of momentum points. From the sufficient number of momenta 
in the Brillouin zone, one can discuss, for example, evolution of the Fermi surface near the 
metal-insulator transition point and momentum dependence of quasiparticle properties in 
detail. Our scheme interpolates high and low energy regions by taking equation of motion 
method combined with the improved DMFT. The high-energy part is systematically taken into 
account by the continued fraction derived from the equation of motion, while the low-energy 
physics is incorporated in the self-consistent scheme of DMFT but here for the self-consistent 
determination of the higher-order self-energy, which now has the momentum dependence in 
contrast to the truncation at the lowest-order self-energy in the original DMFT. 

In order to investigate the fundamental issues we have mentioned, we employ the two 
dimensional Hubbard model with next-nearest-neighbor hopping as is defined by 

w = - Y Y tc ^ c ^ ~Y Y f ' c L c ^ + u Y n ^ n ^- w 

°" (i,j):n.n. & («j'):n.n.n. i 

Here, Ci a (c ia .) is the annihilation (creation) operator of an electron at an atomic site i with a 
spin index a. The number operator on site i, spin a is denoted by ni )(J . The local Coulomb re- 
pulsion is denoted by U. The parameter t and t' are the electron transfers between the nearest 
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neighbor (n.n.) sites and those between the next nearest neighbor (n.n.n.) sites, respectively. 
The summations Ylu j)- n n anc ^ Ylu j)n n n are taken for all the pairs (i,j) between the nearest 
neighbor or between the next nearest neighbor sites, respectively. Here, we set n.n.n. hopping 
parameter t' as t' = — 0.20t. The properties of this model at this value of t' have been inves- 
tigated with several methods. At zero temperature, bandwidth-control (BC) MIT as well as 
filling-control (FC) MIT has been investigated using the path integral renormalization group 
method (PIRG), 13,14 where the first-order metal-insulator transition occurs at a finite value 
of U. At finite temperature, BC MIT has been studied with the correlator projection method 
(CPM), 12 which indicates that the first-order MIT extends to finite temperatures and have a 
critical point at the end of the first-order transition. Finite-temperature study of FC MIT is 
not available and is of great interest. 

In Sec. 2, we briefly summarize the formulation of the correlator projection method. Nu- 
merical results are presented in Sec. 3. Section 4 is devoted to discussions and conclusion. 

2. Formulation 

2.1 Standard Formulation for bandwidth- control problems 

The standard formulation of CPM (for BC problems) is given in previous papers. 12 ' 15 
Here, we just summarize the basic formalism. 

First the equation of motions for operators are solved to give a continued fraction expan- 
sion of the Green's function. We start from a set of operators {cj ct } as the first operator set, 

(2) (3) 

and introduce the second and the third operator set as {(f>)J} and {<p\ a }, respectively. These 
operators are defined through successive operator projection procedure. Here the projection 
of an arbitrary operator O to the n-th operator set {</>' n ^} is defined as 

p n o^Eu°^f n) })( s(n)_1 )z *™' 

i m 

where {•, •} is an anticommutator and (• • •) is the thermodynamic average. The overlap matrix 
for the n-th operator set, ({</>] [ , 4>m }), is denoted by S^ . Then the higher-order operators 
are defined as 

^ + l)_( l _p\ [(j) {n) M 



'I 

with the Hamiltonian of the system TC and commutator [■,•]. Then the Green's function is 
obtained in a continued fraction form as 

G k {iuj n ) = ^ . (2) 

iu n + ii - ei M) - — — ( 2 fc 2) 

iio n + fi - el - Y\ z )[iw n ) 

(1 1) (2 1) (2 2) 

where \i is the chemical potential and the matrices e)' , e,^. , e k ^ for the Hubbard model 
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sf^ = ^(71-^(1 -{n- a }), 



c (2,2) 



[7(1 - (n-a)) - i k - Sfx, 
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(3a) 
(3b) 

(3c) 



where tk is the Fourier transformation of t^, which is defined as 



t 
t' 





(iJ) 

(hj) 
(iJ) 



nearest neighbor 
next nearest neighbor 
otherwise 



We have also introduced a related quantity £/%, which is the Fourier transformation of t^ 
defined as 



1 



,t, 



Uj = (Uj(Si ■ Sj + -mnj - AjAj) ((n- a )(l - (n_ ff ))) , 
where Si and Aj are local spin and pair operators. These are defined as, 

Si = c ia cr a [3Ci[3, 
with Pauli matrices <r, and 

respectively. The equal-time correlation functions of these quantities are calculated with the 
two-particle self-consistent method (TPSC). 16 Finally, 5fi is defined as 

Sfj, = -tij(c\__ a Cj- a (l - 2n i(J )) ((n_ CT )(l - (n_ CT ))) _1 . 

Through these equations, we explicitly included spin indices it's though we hereafter consider 
the spin-symmetric case. Next, in order to evaluate the quantity T^ 2 ', we introduce DMFT 



formulation. Weiss fields for operators c^a and y k 

gf'{iu n ) = 



(2) 



are introduced as 



lU r , 



ln+M -4 1 ' 1 )- e (2,l)g(l)(^ n ) 



(4) 



(5) 



Here, G\ is the local Green's function of the field (j)^ 1 ' , which is formally given below. Then 
E' 2 ) is obtained by a generalized form of the iterated perturbation theory (IPT) of the dy- 
namical mean-field theory 17 as, 






?(°). 



,(0) 



k,k',q c k 



y 1 *,*,,^-^' M»* 



(0) 



,(r), 



(6) 
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where Q k is obtained self-consistently from Eqs. (??),(??) and 

N k iu n + /j,-e k ' ' - £( 2 )(iu; n ) 
Eqs. (4) to (7) constitute a set of self-consistent equations and are solved iteratively until the 
convergence. 

CPM based on the original formulation shown above has achieved a certain success for 
bandwidth-control (BC) problems. 12 For filling-control (FC) problems, however, it shows some 
unphysical results. Namely, the filling changes even when the Fermi level is inside the gap at 
T = 0, while in practice, it should not change within the gap. This failure comes from the 
careless treatment of the chemical potential as we explain in detail below. 

2.2 Reformulation for FC problems 

Now we propose a reformulation devised for FC problems. 

2.2.1 control over the implicit dependence of /i 

In CPM Green's function of the form Eq.(2), its /i dependences not only come from explicit 
ft, but also from implicit /i dependence of £' 2 ) (iio n ) . In order to control /i, all the /i-dependent 
parts have to be controlled consistently. For implicitly /i-dependent function T,( 2 >(iuj n ), this 
is realized by recalculating it with the spectral representation of Green's function. Here, we 
show this procedure for Q, (iuj n ). 

The procedure to determine the chemical potential is the following: 

1. Calculate the spectral function A k (e) of Q k {iu> n ) via Pade approximation. Here, the 
spectral function is defined as A k (e) = — im ^ ( e + * 1 ?)) where n is a small positive number. 

2. Re-calculate the function Q k on the Matsubara axis with desired value of /i by spectral 
representation as 



l^^d 



ioj n + fi — e 

Note that in actual calculation, original Q k to be re-calculated is obtained with a certain 
value of /i|oLD, which is represented as Q k (iuj n ; /i|oLD) • Then the Pade approximation gives 
spectral function with its origin at /u|old as A k (w;/i|oLD) = A k (e = u> + /u|old)- The new 
Q k with chemical potential /i is then obtained as 

Ql \iuj n ;n) = I dwAl (w;/i|oLD)- ; ; • 

J iiv n + fi- /z|old - to 

This Q k is inserted to Eq.(6). 

2. 2. 2 Chemical potential for Weiss fields 

In CPM, Weiss fields for operator Ck a is given as Eq.(4). Here, the chemical potential \i 
appears in the manner just analogous to the lattice Green's function Eq. (2). In fact, however, 
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this form is not obtained by equation of motion analysis like that in the lattice Green's 
function. Therefore it do not have to be necessarily equal to the chemical potential \x in the 
lattice Green's function. Here, in general, Weiss fields are represented with quantity jl as, 

Q k (iwn)fi) = ?YY) ( 8 ) 

iuo n + jl- e\' - e( 2 ' 1 )t/( 1 ) (iuj n ; jl) 

gW(iu n ;P,) = (G^T^iunifi) + S( 2 )(w n ; A))" 1 - (9) 

Here, in Eq.(9), G\ oc {iu> n ] jl) is given by a similar form as the original CPM as 
G[2(iuj n ; fi) 



N ^ iuj n + fl- e( 2 ' 2 ) - E( 2 ) (iio n ; ft) ' 

In the original form of CPM, all jl is set equal to the chemical potential /j, of the lattice 
Green's function. Here, we propose fixing it in a different way in order to satisfy other physical 
requirements. We fix it so as to make the trace of the Weiss field Q k (iu> n ', jl) equal to that of 
the lattice Green's function G k (ito n ; fi), or equivalently, particle number N e ; 



Yl Yl ^ ( iUn '-> ^ = YY1 G k( iuJ n, M) 



k n k n 

The reason for adopting this condition will be given later. 



2.2.3 Properties of the new formulation 

By adding these two reformulation, an important requirement for T = limit is satisfied. 
Requirement: 

If the spectral function is gapped, the result of calculation should not be affected by an 
arbitrary shift of chemical potential within the gap. 

Proof of the statement that the reformulated CPM satisfies this requirement: 
By the first reformulation, chemical potential is completely controlled and the following ob- 
vious relation is satisfied; 

- Y, G k (iv n ; fi) = I deA k {e)f{e - M ). (10) 

P n J 

Here, A k (e) is the spectral function of G k (iuj n ) defined as A k (e) = — ImGfc(e + iO), and func- 
tion f{io) is the usual Fermi distribution function. In case of T — > limit, Fermi distribution 
function becomes a step function. Then the right hand side (RHS) of Eq.(10) becomes 

(RHS) = T deA k (e). (11) 

J-oo 

and it is not affected by the shift of fi to /j,' as long as both of them are within the gap, i.e. 

deA k {e) = / deA k {e) 
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if 

/ deA k {e) = 0. 

Therefore, all the momentum distribution (n&) as well as the total particle number N e is 
unaffected by the shift of chemical potential within the gap. Next, by the second reformulation, 
if the particle number does not change, all Weiss fields are unchanged. This is because fx as 
well as the matrix elements does not change. Note that all the matrix elements such as e^ 1,1 ' 
are only dependent on li and several lattice quantities such as J2j(tijCi a Cj a ) = -^ Sfc*fc( n fc) 
which are unchanged here. Thus the whole results are unaffected. Here, we note that the Weiss 
field is not necessarily gapped when the lattice Green's function has a gap. Therefore the result 
may be changed if one simply set jl equal to li. The idea of fixing the quantity fx different 
from Li to satisfy some physical requirement shares some similarity to the method proposed by 
Kajueter and Kotliar. 18 They proposed a method of applying the iterated pertubation theory 
(IPT) for the dynamical mean-field theory (DMFT) to particle-hole asymmetric models. 19 
Since CPM is an extention of DMFT with IPT, it may also suffer from the problem arising 
from applying IPT in particle-hole asymmetric models. The present prescription solves this 
difficulty. Different solver of CPM could make this procedure unnecessary. 

3. Numerical Results 

In this section, we show the results of our numerical calculations. The following calculations 
have been performed on the 32 x 32 lattice with 1024 (T = 0.05) or 512 (T = 0.10) Matsubara 
frequencies. The spectral functions have been obtained by the analytic continuation using 
Pade approximation. Throughout this paper the energy unit is set as t = 1. 

3.1 Phase diagram 

First, we show phase diagram of the Hubbard model at fixed temperatures. At T = 0.05 
and 0.10, the li-U phase diagram is obtained as Fig. 1 (For an enlarged figure see also Fig. 2) . 
Here, first-order transitions appear with a certain range of coexsistence of two solutions. The 
phase transition obtained here contains Lifshitz transitions, which is the transition caused by 
the topology change of the Fermi surface. In Fig.l, roughly speaking, the equilibrium solution 
shows a first-order transition from metals with a large Fermi surface to insulators. However, 
if one carefully examines the metastable solutions obtained from the adiabatic continuation 
of the insulating phase, it turns out that the Fermi surface changes from large to small or 
vice versa as is shown in Fig. 3. Therefore, it turns out that the Lifshitz transition occurs 
simultaneously with the metal-insulator transition. The small Fermi surface is preempted in 
the equilibrium because of the first-order transition. We refer to the phase with the small 
Fermi surface as S-phase and that with the large Fermi surface as L-phase. We also use the 
expression such as S-metallic phase and L-metallic phase in order to describe metallic phase 
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Fig. 1. Phase diagram of the 2D Hubbard model at t' = —0.20 in the plane of \x and U at T — 0.05 
(upper panel) and at T = 0.1 (lower panel). Solid lines indicate the simultaneous first-order Lifshitz 
transition and metal- insulator transition (MIT), while the dashed lines indicate the parameters 
with which the metastable solution becomes unstable. Metastable insulating solution exist in the 
regions sandwiched by the solid and the dashed lines. The squares are the critical points of the 
Lifshitz transitions, or equivalcntly, the cndpoints of the first-order Lifshitz transitions. Note that 
the dash-dotted lines are drawn between the square and the circle in small segments of first-order 
Lifshitz-transition line and indicate first-order Lifshitz transitions without MIT. See Fig. 2 for more 
detail. Circles, which separate solid lines and the dash-dotted lines are also explained in Fig. 2. 
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Fig. 2. Detailed view of the [i-U phase diagram at T = 0.10, enlarged from Fig. 1. Here, as the 
same as Fig. 1, the solid and dash-dotted lines indicate the first-order Lifshitz transition, where 
the solid line accompanies MIT while the dash-dotted line does not. Thick dotted line indicates 
metal-insulator crossover lines. On this line, the bottom of the upper Hubbard band (UHB) crosses 
the Fermi level and the system crosses over between metallic phase and insulating phase. Note 
that such crossing of QP dispersion does not indicate a transition since the system is at a finite 
temperature. On the other hand, along the solid lines above the circles, the crossing occur dis- 
continuously and one can define metal-insulator transitions. The thin dashed lines indicate the 
parameters with which the metastable insulator solution becomes unstable. 



with small Fermi suface and that with large Fermi surface, respectively. Although the S-phase 
in Fig. 3 stays as a metastable one, we will show below that the S-phase indeed becomes the 
equilibrium in the region between the square and the filled circle in Fig.l. 

We note that Lifshitz transitions as well as metal-insulator transitions are originally tran- 
sitions in the zero-temperature limit. However, they are extended to finite temperature if these 
transitions are of the first order. Finite temperature Lifshitz (metal-insulator) transitions are 
defined as the first-order transitions from phase A to phase B where the zero-temperature 
extention of the phase A and the phase B are separated by the Lifshitz (metal-insulator) tran- 
sitions of the original definitions. The spectral density distribution changes discontinuously 
at first-order Lifshitz transitions. Therefore they accompany metal-insulator transitions if the 
spectral density changes from metallic to insulating. In the phase diagram in Fig. 1 as well as 
in Fig. 2, solid lines above the circles show such kind of transitions. 
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We consider the region above the solid circle in Fig.l in more detail, where the first- 
order Lifshitz transition occurs simultaneously with MIT. Typical calculated results for the 
density and the grand canonical potential at a fixed value of U, U = 3.50 and sweeping 
/j, in different directions are shown in Figs. 4 and 5. Here, a clear hysteresis is observed. 
One series of solutions is obtained by sweeping n from the lowest values of fj, to n = U/2, 
while the other series of solutions is obtained by sweeping fi in the opposite direction. As 
will be explained below, each corresponds to the sweep from the hole-doped (electron-doped) 
metallic phase to the insulating phase and vice versa, respectively. To clarify the nature of 
each solution, we take the electron-doping side for example, and plotted the Fermi surface of 
each solution in Fig. 3. Here, each Fermi surface is obtained as the zero-energy section of the 
QP dispersion. The solutions obtained through sweeping [i in the increasing order, which are 
metastable solutions as we will explain below, are solutions with the small Fermi surfaces ( 
S-phase ). Those obtained through sweeping [i in the decreasing order, which are equilibrium 
solutions, are those with larger Fermi surfaces ( L-phase ). We compare the thermodynamic 
potential S1|gce = F — /j,(N e ), where F is the Helmholtz's free energy, for the both solutions. 
As shown in Fig. 5, the L-phase has lower potential compared to the S-phase in most part of 
the coexistence region, and the equilibrium phase boundary is very close to the point where 
the L-phase solution becomes unstable. 

In this parameter region, spectral density of the metastable region of the S-phase changes 
from the metallic to insulating as the system approaches half filling. It is insulating at the 
point where the true phase transition occurs. Therefore, the transition here accompanies a 
metal-insulator transition. 

Now we discuss the detailed view of the [i-U phase diagram shown in Fig. 2 in detail. 
Along the dash-dotted lines below the circles, first-order Lifshitz transition occurs from S- 
metallic to L-metallic and the metal-insulator transition becomes just a continuous crossover 
separated from the Lifshitz transition. Here, we investigate the behavior around the critical 
point at fixed U = U c for several choices of ji. As is shown in Fig. 2, as the chemical potential 
increases from the insulating region at U = U c , the system crosses over from the insulating 
phase to the the S-metallic phase. Then it crosses the critical point, and finally, enters the 
L-metallic phase. We have shown the evolution of the Fermi surface in Fig. 6. Within the 
rigid band picture, the insulating phase crosses over to the S-metallic phase as the chemical 
potential jjl crosses the bottoms of the upper Hubbard band (UHB), which are the M-points 
(here defined as the momentum (ir, 0) and its equivalent points). In the S-metallic phase, 
the section of the Fermi surface is disconnected pockets around the M points. It changes to 
connected large one at the Lifshitz transition point as the Fermi level crosses the X-points 
(here defined as the momentum (tt/2, 7t/2) and its equivalent points). In this intermediate 
region, the S-metallic phase is indeed the equilibrium. Apparently, the S-metallic phase does 
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Fig. 3. Fermi surfaces of solutions at fixed U = 3.50 and different values of /j. — J7/2. Starting from 
(f), in the sweep of decreasing /i, we obtain (f)-(d)-(b)-(a) in a sequence, while in the sweep of 
increasing /i, we obtain (a)-(c)-(e)-(f). Here, at [i — U/2 = 0.85 and at [i — U/2 = 1.05 of the /it- 
increasing sweep, the solutions are insulating and have no Fermi surface. There is a discontinuous 
transition between (b) and (a) in the former sweep, while between (e) and (f) in the latter sweep. 
The transition (c) to (f) is not a metal-insulator transition, but a Lifshitz transition from the 
small Fermi surface (S-phase) to the large connected Fermi surface (L-phase). At the transition 
(b) to (a), the Lifshitz transition accompany a metal-insulator transition. As is mentioned in the 
text, the solutions obtained through decreasing /i are the true equilibrium 

not satisfy the Luttinger sum rule. 20 Since the Mott gap is formed without any long-ranged 
order in this formalism, the Luttinger volume cannot be satisfied near half filling unless this 
S-metallic phase is preempted by a strong first-order transition. The violation of Luttinger 
sum rule in underdoped region was also suggested in cluster DMFT. 21 At the transition in the 
electron-doped side, appearance of Fermi pockets in the underdoped region similar to those 
in the S-metallic phase was shown in certain approximations. 22 ' 23 In their cases, however, the 
symmetry breaking by the antiferromagnetic order was assumed as the starting point, while 
we have not assumed the symmetry breaking in the present study. At nonzero temperatures, 
the S-metallic phase is in any case adiabatically connected with the insulating phase and a 
strict statement on the Luttinger sum rule is neither possible nor meaningful. Although it is 
remarkable that the apparent small Fermi pocket is detected in our study at low but finite 
temperatures, it does not necessarily exclude that Lifshitz transition and MIT always merge 
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Fig. 4. Plot of filling (n) at T = O.fO, U = 3.50 as a function of the chemical potential \l. The solid 
line with /i — U/2 < (fj, — U/2 > 0) was obtained by sweeping /j, from the lowest (highest) values 
of /i to (i — U/2 = 0, while the dashed line was obtained by sweeping /i in the opposite direction. 
Note that the full line corresponds to the sweep from the metallic to the insulating phase, and the 
dashed line, opposite. 
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Fig. 5. Plot of the grand canonical potential f^GCE for the solutions shown in Fig. 4. Solid and dashed 
lines show the sweep in the different directions as is shown in Fig. 4. The result shows that the 
solutions along the solid line are the equilibrium ones within error bars. 

at strictly zero temperature, thereby the S-phase is preempted by the first-order transition. 
However, the critical temperatures of Lifshitz and metal-insulator transitions may be lowered 
to zero by enhancing quantum fluctuations, where the first-order transitions disappear. If the 
both transitions become continuous at zero temperature, we do not know how to make them 
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compatible with the Luttinger sum rule, if the quasiparticle weight remains nonzero through 
the transitions as we clarify later. The gap formation without the symmetry breaking shares a 
similarity to the renormalization group results 24 The region of this separation between metal- 
insulator and Lifshitz transitions may be more extended for larger t' with larger region of 
S-metallic phase. 

The S-metallic phase sandwiched by the Lifshitz transition and the metal-insulator transi- 
tion in the underdoped region may show unusual metallic behavior far from the simple Fermi 
liquid. Near the S-metallic phase, the small Fermi volume as in Fig. 6(c) has very anisotropic 
Fermi surface and the inner circle has faint weight as compared to the outer surface. When 
the faint inner surface becomes not observable in the experimental resolution, this may lead 
to a large Fermi volume obtained from the outer curve. The volume quickly approaches that 
expected from the Luttinger theorem with evolving doping. At the same time, this anisotropy 
in the case of Fig. 6(b) may lead to observation of an arc structure of the Fermi surface only 
consisting of a part of the Fermi surface. In the hole doped case, this makes the arc structure 
around (7r/2, tt/2) as is observed in angle resolved photoemission experiments of the cuprate 
superconductors. 6 ' 7 The Luttinger sum rule may be approximately satisfied for relatively small 
t', while it may be strongly violated in the S-phase as well as in the situation like Fig. 6(c) 
with increasing t' at a relatively small U region. 

3.2 Determination of critical points 

As is shown in Fig. 1 there are critical points at the ends of the first-order phase transition 
lines. The location of the critical points are carefully determined by an extrapolation procedure 
which is described below : We have obtained the first-order phase transition lines. Beyond 
the Lifshitz critical points, which are shown as squares in Fig. 1, there are crossover lines, 
though the crossover lines are not shown explicitly in the figure. The crossover lines are the 
extensions of the first-order phase transition lines and defined as the maximum of the charge 
susceptibility Xc = q 1> - The precise location of the critical points are found along this line. 
From the first-order transition side, it is found by extrapolating the first-order jump to zero. 
From the crossover side, it is found by extrapolating the inverse of the charge susceptibility 
to zero. An example of this extrapolation is shown in Fig. 7. 

3.3 Analysis of the critical points 

One of the major subjects of this research is to investigate the nature of these critical 
points. Now we analyze them in detail. The critical points are the points where the charge 
susceptibility \c diverges. The charge susceptibility is related to the spectral density function 
Ak{oj) as 

d(n) d f 2 



Xc = -^f = Q^ I d 2 k J deA k (e)f(e-v,) (12) 
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Fig. 6. The Fermi surface around the critical point at fixed value of U; U = U c = 1.05, and different 
values of \x\ \x - U/2 = 0.1085, 0.1485, 0.1585, and 0.1685, respectively. Note that at /j, - U/2 = 
0.1085, Fermi surface does not exist since it is in the insulating phase. As the chemical potential 
increases from /i — U/2 = 0.1085, the system crosses over from insulating phase (fi — U/2 = 0.1085) 
to metallic phase with small Fermi surfaces ( S-metallic ) (fi — U/2 = 0.1485). It then crosses the 
critical point (p, — U/2 = 0.1585), and finally metallic phase with large connected Fermi surface 
( L-metallic ) (/j, — U/2 = 0.1685). The weight of the inner circle in (c) and (d) is much smaller 
than that of outer open curves. With further evolution of doping beyond (d), the inner circle in 
(d) fades out and disappears which leaves the single Fermi surface evolved from the outer curves. 



Now we introduce the concept of quasiparticle (QP). We label each QP by its momentum k 
and index a, where the latter is introduced in order to distinguish different QP peaks with 
the same momentum k. The QP with label k, a is characterized by its pole energy eta-, an d 
the pole weight Z^ a . Then the spectral function A^^co) is devided into QP contribution and 
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Fig. 7. Singularities around the electron-like critical point. Data points are taken along the first- 
order Lifshitz-transition line and its extension, crossover line (see text). From the continuous 
transition side, from lower U, minimum values of inverse charge susceptibility (dn/d^,)^ 1 are 
plotted with circles, while from the discontinuous transition side, jumps of the filling An at the 
Lifshitz transition points are plotted with triangles. Here, the circles are extrapolated linearly to 
the zero-point at U — 1.05. The triangles do not show linear behavior, and are monotonic, which 
also decreases to zero near U c — 1.05. These suggest the exponents 7 = 1 and /? < 1, which are 
consistent with the mean-field Ising values 7=1 and j3 = 1/2, although quantitative estimates 
are left for future studies. Note that the exponents are defined by (dn/d^,)^ 1 oc \U — U c \' y and 
Anoc \U -U c f. 



the remaining incoherent part as 

A fc (e) = J2 Z ^(e ~ e ka ) + ^(e). (13) 

a 

After inserting Eq.(13) to Eq.(12) we obtain the following expression for Xc ', 
Xc=^jd 2 kY J Z ka f(e ka -n) + ^jd 2 kj de4 ncoh (e)/(e - M ) 

= J « t„ * E ^i|(0 + ± / A / *«/(e - „>. (U) 

If the concept of QP plays an important role, the first term in the last line of Eq. (14) becomes 
dominant, and the divergence of the charge susceptibility occurs through the term 



!Ve fco 



i.e., 



flattening of the QP dispersion. 

Next we show the procedure which we have taken in order to derive the above quantities 
from the numerical calculation. 
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1. Quasiparticle (QP) dispersion e ka 

We calculate spectral functions with small imaginary part r\ as A2(oj; /i) = — ImGk(u>+ir]; //). 
QP peaks are denned as the peak structure in A^(uj; fi) with intensity larger than 1% of its 
maximum value — . 

2. Quasiparticle weights 

We define the QP pole weight Zqp(A;, a) as an integration of spectral function around the QP 

pole; 

/•A/2 

ZQ P (k,a)= / dJA k (J + e ka ). 
J -A/2 

Here, A is the width of energy window. We fix it as A = T. This quasiparticle weight is 

computed for the pole, which can be adiabatically connected to the quasiparticle at the Fermi 

level in the metallic side. A related quantity is their Fermi- level contribution Zps(k), which 

is defined as an integration of spectral function around the Fermi level; 

/•A/2 
Z FS (k) = / dw'A k (u'). 

J-A/2 

3. Effective Fermi momentum 5k c p(5ii) 

Here, we define the Fermi surface by the zero energy section of the QP energy dispersion 
surface. Around the critical point of the Lifshitz transition on the L-phase side, we introduce 
the following quantities to estimate the chemical potential dependence of the evolution of the 
Fermi surface: At the critical point, the Fermi surface has a topological singlularity. We define 
the momentum of this point as k-p c . I n the L-phase side, the Fermi surface expands around 
&Fc- We measure the Fermi momentum at chemical potential fi, &f(aO from the critical k Fc - 
This quantity is a function of chemical potential measured from the critical chemical potential 
fi c . Now we define 5/j, = \fi — fi c \, and 5k c p(5fj>) = kp(^) — kp c . We note that in the rigid band 
limit, 5k c p(5ii) is equivalent to the QP dispersion at the critical point, whereas it is not in 
case the QP dispersion changes with the chemical potential. 

Through the analysis, xy-symmetry is always conserved. Therefore, as for momentum 
space quantities, xy-symmetric points are equivalent. Hereafter, we refer to (ir, 0) and its 
equivalent points in the momentum space as M-points, and (tt/2, tt/2) and its equivalent 
points in the momentum space as X-points. With these preparations, we analyze one of the 
obtained critical points. We choose the critical point at T = 0.10 and study transitions to 
electron-like Fermi surface. 

3.4 QP dispersion around the critical point 

In the critical endpoint of the Lifshitz transition at finite temperature, the compressibility 
diverges as we see above. To get insight into the origin of this divergence of the charge suscep- 
tibility Xc at the Lifshitz critical point, we examine the quasiparticle dispersions. The QP dis- 
persion around the critical point is plotted in Fig. 8. It shows a saddle point at (7ir/16, 7ir/16), 
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Fig. 8. QP dispersion of the UHB at U = U c for several choices of \x; \x — U/2 — 0.1485, 0.1585, and 
0.1685, respectively. The QP pole energy of each k-point is shown with color. The white region 
corresponds to fc-points with no sharp QP pole in UHB. Note that blue, purple, and gray regions 
correspond to Cf- < 0. At [i — U/2 = 0.1585, the critical point, X' points are at Fermi level, tf. = 
and are the topological critical points that connect two hole pockets formed around M-points. 
X'-points are also saddle points of QP dispersion and relatively flat dispersion is found near them. 



which leads to a flat dispersion near the X-point. We refer to this saddle point and its equiv- 
alent points as X' points. The QP dispersion at critical value of /j, around this X' point is 
quadratic in (1, 1) direction. Thus the divergence of the charge susceptibility is explained 
within the rigid band picture as the logarithmic divergence at the saddle point; 

d{n) If ..1 

" N 



<9/x 



dk- 
e fc =o |ve fcc 



(15) 



In order to further investigate the flattness of the QP-dispersion around the X' points, we 
plotted the quantity 5k c p(5fi) around the X' points. As we described above, this quantity 
shows the chemical potential dependence of the evolution of the Fermi surface and expected 
to be more directly related to the divergence of the charge susceptibility. We have derived 
/c-exponent z as {5k e p) z = 5fi as is illustrated in Fig. 9. Here the obtained z is z ~ 2.5 > 2, 
which also explain the divergence of the compressibility by the saddle point singularity of 
QP dispersion. Note that we have fitted the dispersion from —(1, 1) direction data only since 
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in that direction, 5kp(5/j,) extends more. We have also discarded data for Skp < y 2, since 
it depends on the interpolating procedure. 25 Although we need much larger system size to 
determine the precise exponent, z larger than 2 obtained within the present accuracy supports 
the enhancement of flattening of the effective band dispersion beyond the simple van Hove 
singularity. 

3.5 QP weight around the transition 

One of the central issues in this paper is to clarify the nature of the metal-insulator tran- 
sition and fate of the quasiparticles at the transition. In §1, we have introduced two different 
pictures of MIT. In order to clarify the origin and criticality of MIT, we have calculated 
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pole weights of quasiparticles around the critical Lifshitz points as well as around the metal- 
insulator transition. As is shown in Fig. 10, QP weight at the Fermi level decreases to zero 
as the system moves from the metallic to the insulating phase. Here, we have plotted the 
change of the QP weights across the transition line at three points with different value of 
U; U = 1.05, U = 1.30, and U = 3.50, which are the electron-like critical point of Lifshitz 
transition, the end of MIT point, and a point in the strong-coupling region, respectively. The 
Lifshitz transition is continuous at the first point and of the first order at the latter two. Note 
that in the plot of U = 1.30, the energy of the lowest energy QP in the UHB becomes zero. QP 
weight at the Fermi surface remains nonzero even in slightly insulating side of the transition 
due to finite temperature effects. 

Clearly, the QP weight at Fermi level vanishes as the system changes to insulating side. 
On the other hand, the QP weights of each QP poles remain finite through the transition. 
Although the weight jumps at the first-order metal-insulator transition points, the QP weight 
Zqp remains nonzero in the insulating side even through the first-order metal-insulator tran- 
sition. The QP weights Zqp at points other than M point show similar behavior. 

The result indicates that the metal-insulator transition is driven by the shift of Fermi sur- 
face out of the QP bands ( LHB or UHB ) rather than vanishing of QP weights irrespective of 
the order of MIT and irrespective of the involvement of the Lifshitz transition. It supports the 
quasiparticle rigid-band picture of MIT rather than the quasiparticle- weight vanishig picture 
of MIT. QP excitations are still well-defined around the transition although their dispersions 
change in a non-trivial manner around the transition. Analysis of QP dispersions is crucial for 
the study of Mott-criticality. This result is consistent with quantum Monte Carlo calculation 
of a single hole in the t-J model. 26 These results indicate the validity of recently reported QP- 
based analysis of the Mott criticality. 27-29 It is remarkable that the higher-order correction of 
DMFT to include the momentum dependence presented here as CPM seriously modifies the 
character of the Mott transition from the quasiparticle- weight vanishing type in the original 
DMFT to the quasiparticle rigid band type. 

4. Discussion and Conclusion 

Using CPM, which takes into account momentum dependence of the self-energy beyond 
the dynamical mean-field theory, we find the Lifshitz transition playing a substantial role near 
the MIT. At the Lifshitz transition, the Fermi level crosses the singular point in quasiparticle 
dispersion. This causes a large shift in the QP dispersion and, within a certain parameter 
region, yields the first-order MIT. Here, the Lifshitz transition itself causes a large shift in 
the momentum distribution. However, the self-consistent scheme of CPM relates this shift to 
the shift in the QP dispersion. Interestingly, such transition of Fermi surface toplogy is also 
obtained on the same model, with several different calculation methods. Recent study with 
the cellular DMFT method reports the change of the Fermi surface topology near half-filling 
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Fig. 10. Plot of QP weights as a function of chemical potential /i. Left-top panel, around the electron- 
like critical point ( U = 1.05 ), right-top panel, around the MIT point (U = 1.30), bottom panel, 
strong coupling region with discontinuous Lifshitz and metal- insulator transition line ( U = 3.50 ). 
In each panel, largest value of Zps out of the momentum space, .Zfsmax, as well as Zqp and Zps 
for representative momenta, namely at M point, are plotted against chemical potential /i. Note that 
QP at the M-points is the lowest energy pole of the UHB and play key role for the metal-insulator 
transition or crossover. Here, the quantity ^fsmax shows small non-monotonic behavior, though 
we attribute this to an artifact by finite-size effects coming from discrete momentum points. The 
overall feature does not suffer from the size effects. In all the cases, the quasiparticle weight at the 
M point remains nonzero through the transition. 



although the relation to MIT has not been studied. 30 

However, we also note that all the above results including the present one stays ulti- 
mately at the mean-field level of approximations. By considering the quantum fluctuations 
more precisely, the first-order transition and the critical temperature might be suppressed. 
The compressibility divergence then may occur only at zero temperature in a wide region of 
the metal-insulator boundary. This is nothing but the marginally quantum critical behavior 
recently discussed by one of the authors. 2 This possibility is supported in the path-integral 
renormalization group study, 14 where the quantum fluctuations are fully taken into account. 

Clear Lifshitz transitions are also observed in a recent study with the Hartree-Fock ap- 
proximation (HFA). 31 In case of HFA, however, it occurs in much more heavily doped (either 
electron-doped or hole-doped) regime and is rather independent of the metal-insulator tran- 
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sition. If the Luttinger sum rule 20 holds, the critical shape of the Fermi surface determines 
the critical values of filling. When the momentum of the van Hove singularity is far from the 
folded Broulline-zone boundary for the insulator, the doping concentration required to cause 
the Lifshitz transition becomes large. This is indeed the case of the Hartree-Fock approxima- 
tion at large t'/t. The doping concentration xlfst at the Lifshitz transition is given by the 
area of the Fermi surface pocket, since the pocket vanishes at the metal-insulator transition. 
On the other hand, in CPM, the hole/electron pockets consist of quasiparticles, with a small 
factor Zqp. Therefore the filling at the Lifshitz transition point may become very close to the 
filling at the MIT (pockets-vanishing) point. They in fact coincide and xlfst becomes zero 
for larger U/t. This may also be interpreted by the renormalization of the Fermi surface shape 
to the perfectly nested one by the self-energy effect. 

When we consider the long-range Coulomb repulsion in addition to the Hubbard model, 
the first-order transition is suppressed even when the quantum fluctuations are not seriouly 
considered. In this situation, the region of the first-order jump in the density (namely, the 
region of lightly doped Mott insulators at large U) becomes a nontrivial phase, where the 
uniform electron density is unstable but the complete phase separation is suppressed. Quantum 
fluctuations may again smear out this peculiar phase and lead to a crossover with an unusual 
metallic behavior. 

In conclusion, we have proposed a revised form of CPM to investigate filling-control metal- 
insulator transition. This formalism of CPM allows incorporating momentum dependence of 
the self-energy beyond the dynamical mean-field theory with large number of momenta re- 
solved. This makes it possible to study evolution of quasiparticle spectra as well as phase 
diagram near the Mott transition in detail. We have studied the filling-control metal-insulator 
transition in the two-dimensional Hubbard model. Within CPM, we conclude that the quasi- 
particle survives on the edge of the Mott gap through the metal-insulator transition. The 
Mott transition is characterized by the Fermi level being separated from the quasiparticle 
dispersion and entering the gap by keeping nonzero quasiparticle weight. In this sense, the 
rigid band picture is ultimately valid on the verge of the transition, where the preformed Mott 
gap and quasiparticle are both retained through the transition. At the same time, the actual 
behavior of the QP dispersion is substantially modified from the original simple rigid band 
picture. In a part of the phase diagram, the metal-insulator transition occurs simultaneously 
with the Lifshitz transition as a combined first-order one. On the other hand, it also shows 
a region where the two transitions are separated, between which an unusual metallic phase 
with violation of the Luttinger sum rule appears to be stabilized. 
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